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Abstract 

In contrast to the conventional wisdom that scale-free networks are prone to epi- 
demic propagation, in the paper we present that disease spreading is inhibited in 
fractal scale-free networks. We first propose a novel network model and show that 
it simultaneously has the following rich topological properties: scale-free degree dis- 
tribution, tunable clustering coefficient, "large-world" behavior, and fractal scaling. 
Existing network models do not display these characteristics. Then, we investigate 
the susceptible- infected-removed (SIR) model of the propagation of diseases in our 
fractal scale-free networks by mapping it to bond percolation process. We find an 
existence of nonzero tunable epidemic thresholds by making use of the renormal- 
ization group technique, which implies that power-law degree distribution does not 
suffice to characterize the epidemic dynamics on top of scale- free networks. We argue 
that the epidemic dynamics are determined by the topological properties, especially 
the fractality and its accompanying "large-world" behavior. 
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1 Introduction 



In recent years, there has been much interest in the study of the structure and 
dynamics of complex networks [T|2||3|H] . One aspect that has received consid- 
erable attention is the epidemic spreading taking place on top of networks [5], 
which is relevant to computer virus diffusion, information and rumor spread- 
ing, and so on. In the study of epidemic spreading, the notion of thresholds is 
a crucial problem since it finds an intermediate practical application in disease 
eradication and vaccination programs [SfT] . In homogeneous networks, there is 
an existence of nonzero infection threshold, if the spreading rate is above the 
threshold, the infection spreads and becomes endemic, otherwise the infection 
dies outs quickly. However, recent studies demonstrate that the threshold is 
absent in heterogeneous scale- free networks [8115111 Of 11] . Thus, it is important 
to identify what characteristics of network structure determine the presence 
or not of epidemic thresholds. 

To date the influences of most structural properties on disease dynamics have 
been studied, which include degree distribution [9|10|lllj . clustering coeffi- 
cient [T2], and degree correlations [T3] . However, these features do not suffice 
to characterize the architecture of a network [TJ]. Very recently, by introduc- 
ing and applying box-covering (renormalization) technique, Song, Havlin and 
Makse found the presence of fractal scaling in a variety of real networks [T5][T6] . 
Examples of fractal networks include the WWW, actor collaboration network, 
metabolic network, and yeast protein interaction network (20] • The fractal 
topology is often characterized through two quantities: fractal dimension ds 
and degree exponent of the boxes dk, both of which can be calculated by the 
box-counting algorithm [TTfT8] . The scaling of the minimum possible number 
of boxes Nb of linear size £b required to cover the network defines the fractal 
dimension ds, namely Nb ~ £~b b - Similarly, the degree exponent of the boxes 
dk can be found via k B (£B)/kh u b ~ ^ k i where &b(£b) is the degree of a box 
in the renormalized network, and khub the degree of the most-connected node 
inside the corresponding box. Interestingly, for fractal scale-free networks with 
degree distribution P(k) ~ A; -7 , the two exponents, ds and d k , are related to 
each other through the following universal relation: 7 = 1 + d B /d k [15] . 

Fractality is now acknowledged as a fundamental property of a complex net- 
work [H]. It relates to a lot of aspects of network structure and dynamics 
running on the network. For example, in fractal networks the correlation be- 
tween degree and betweenness centrality of nodes is much weaker than that 
in non-fractal networks [TU]. In addition, several studies uncovered that frac- 
tal networks are not assortative [T61I20II21] . The peculiar structural nature of 
fractal networks make them exhibit distinct dynamics. It is known that frac- 
tal scale-free networks are more robust than non-fractal ones against mali- 
cious attacks on hub nodes [ToTTST] . On the other hand, fractal networks and 



2 



their non-fractal counterparts also display disparate phenomena of other dy- 
namics, such as cooperation [23.22J , synchronization [21J, transport |24|, and 
first-passage time [2"t3f2~o] . Despite of the ubiquity of fractal feature and its 
important impacts on dynamical processes, the dynamics of disease outbreaks 
in fractal networks has been far less investigated. 

In this current paper, we focus on the effects of fractality on the dynamics 
of disease in fractal scale-free networks. Firstly, we propose an algorithm to 
create a class of fractal scale-free graphs by introducing a control parame- 
ter q. Secondly, we give in detail a scrutiny of the network architecture. The 
analysis results show that this class of networks have unique topologies. They 
are simultaneously scale-free, fractal, 'large- world', and have tunable cluster- 
ing coefficient. Thirdly, we study a paradigmatic epidemiological model [61T7] . 
namely the susceptible-infected-removed (SIR) model on the proposed fractal 
graphs. By mapping the SIR model to a bond percolation problem and using 
the renormalizat ion-group theory, we find the existence of non-zero epidemic 
thresholds as a function of q. We also provide an explanation for our findings. 



2 Network construction and topological properties 

This section is devoted to the construction and the relevant structural proper- 
ties of the networks under consideration, such as degree distribution, clustering 
coefficient, average path length (APL), and fractality. 



Fig. 1. (Color online) Iterative construction method of the fractal networks. Each 
iterative link is replaced by a connected cluster on the right-hand side of the arrow. 
The red link is a noniterated link. 

2.1 Construction algorithm 

The proposed fractal networks have two categories of bonds (links or edges): 
iterative bonds and noniterated bonds which are depicted as solid and dashed 
lines, respectively. The networks are constructed in an iterative way as shown 
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in Fig. [H Let F t (t > 0) denote the networks after t iterations. Then the 
networks are built in the following way: For t = 0, F is two nodes (vertices) 
connected by an iterative edge. For t > 1, F t is obtained from F t _i. We replace 
each existing iterative bond in F t _i either by a connected cluster of links on 
the top middle of Fig. [1] with probability q, or by the connected cluster on the 
bottom right with complementary probability 1 — q. The growing process is 
repeated t times, with the fractal graphs obtained in the limit t —>■ oo. Figure [2] 
shows a network after three-generation growth for a specific case of q = 0.6. 




Fig. 2. (Color online) Sketch of a network after three iterations for the particular 
case of q = 0.6. 

Next we compute the numbers of total nodes and edges in F t . Let L v (t), 
Li(t) and L n (t) be the numbers of nodes, iterative edges, and noniterated 
edges created at step t, respectively. Note that each of the existing iterative 
edges yields two nodes and four new iterative edges; at the same time this 
original iterative edge itself is deleted, which means that Li{t) is also the total 
number of iterative edges at time t. Then we have L v (t) = 2 Li(t — 1) and 
Li(t) = 4Li(t — 1) for all t > 0. Considering the initial condition L«(0) = 1, 
one can obtain Li(t) = 4* and L v (t) = 2 x 4* -1 . Thus the number of total 
nodes N t present at step t is 

- 2 x 4* -I- 4 

On the other hand, at each construction step, each of the existing iterative 
edges may yield one noniterated link with probability 1 — q, so the expected 
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value of L n (t) is (1-q) L^t-l) fort > 1, i.e., L n (t) = (1-q) 4* _1 . Therefore, 
the total number of edges E t present at step t is 

g , = L, w + E^) = (4 " 9)4 '" (1 " 9) . 

The average node degree after £ iterations is (k) t 
4 — g in the infinite t limit. 

2.2 Degree distribution 

When a new node u enters the system at step t u (t u > 1), it has two iterative 
edges. At the same time, with probability 1 — q one noniterated edge is created 
and linked to node u. Let Li(u,t) be the number of iterative links emanated 
from node u at step t, then Li{u,t u ) = 2. Notice that at any subsequent step 
each iterative edge of u is broken and generates two new iterative edges linked 
to u. Thus Li(u, t) = 2 Li(u, t - 1) = 2*"*" +1 . 

We define k u (t) as the degree of node u at time t, then we have 

k u (t) = 2 t ~ tu+1 , with probability q, 

k u (t) = 2 t ~ tu+1 + 1, with probability 1 — q, 

where the last term 1 in the second formula represents the noniterated link 
connected to node u. For the initial two nodes created at step 0, neither of 
them has a noniterated link, both nodes have a degree of 2*. 

Equation ([3]) indicates that the degree spectrum of the networks is not con- 
tinuous. It follows that the cumulative degree distribution [2"T1I2"%] is given by 
-Pcum(^) — where N t ^ k is the number of nodes whose degree is not less than 
k. When t is large enough, we find P cum (^) ~ k~ 2 . So the degree distribution 
follows a power-law form with the exponent 7 = 3, which is independent of q. 
The same degree exponent has been obtained in the famous Barabasi- Albert 
(BA) model 0. 

2.3 Clustering Coefficient 

The clustering coefficient [29] of a node u with degree k u is the probability 
that a pair of neighbors of u are themselves connected, which is given by 
C(k u ) = 2b u /[k u (k u — 1)], where b u is the number of existing connections 



(2) 

= which approaches 
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between the k u neighbors of u. For our networks, the clustering coefficient 
C(k) for a single node with degree k can be evaluated exactly. Note that 
except for those nodes born at step t, all existing links among the neighbors of 
a given node are noniterated ones, whose number is ease to calculate. For the 
initial two nodes, the expected existing noniterated links among the neighbors 
is (1 — q)2 t ~ 1 . For each of those nodes created at step (0 < (p < t), there 
are average (1 — q)2 t ^^ > noniterated links among its neighbors. Finally, for 
the nodes generated at step t, some of them have a degree of 3, the number 
of links between the neighbors of each is 2; the others have a degree of 2, 
their clustering coefficient is zero. Thus, there is a one-to-one correspondence 
between the clustering coefficient C(k) of the node and its degree k: 



C{k) 



(l-g)/(Jfe-l) for k = 2 m (2 <m<t), 

(l-q)/k for k = 2 m + 1(2 < m < t), 

for k = 2, 

2/k for k = 3. 



(4) 



Using Eq. (J3J), we can obtain the clustering C t of whole the network at step 




Fig. 3. (Color online) The average clustering coefficient C of the whole network as 
a function of q. It can be tunable systematically by changing q. The squares are 
the simulation results, and the line represents the analytical expression given by 
Eq. ©. 

t, which is defined as the average clustering coefficient of all individual nodes. 
Then we have 



L v (0) ^ qL v (r) 
2* - 1 r ti K r - 1 



+ g(l-g)L„(r) , 2L v (t) 

r=l 



K r + 1 



K t + 1 



(5) 
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where K r = 2* _r+1 and K t = 2. In the infinite network order limit (N t — > oo), 
C{ converges to a nonzero value C as a function of parameter q. In the two 
extreme cases of q = and q = 1, C are and 0.5435 [50] . respectively. 
When 5 increases from to 1, C grows from to 0.5435, see Fig. [3J Thus, 
the parameter q in our model introduces the clustering effect by allowing the 
formation of triangles. Furthermore, the relation between C and q is almost 
linear, as depicted in Fig. [31 

2.4 Average path length 

Shortest paths play an important role both in the transport and communica- 
tion within a network and in the characterization of the internal structure of 
the network [31f 32 33 34]. Let dij represent the shortest path length from node 
i to j, then the average path length (APL) d t of F t is defined as the mean of 
d^ over all couples of nodes in the network, and the maximum value D t of d^ 
is called the diameter of the network. 

For general q, it is difficult to derive a closed formula for the APL d t of 
network F t . But for two limiting cases of q = and 5 = 1, both the networks 
are deterministic ones, we can obtain the analytic solutions for APL, denoted 
as g?° and dj, respectively. In the large t limit, cf° ps |^2* [30] and d] rs |j2* [23] . 
On the other hand, for large t limit, N t ~ 4', so both d® and d\ grow as a 
square power of the number of network nodes. 

By construction, it is obvious that d® and d\ are the lower and upper bound of 
APL for network F t , respectively. As a matter of fact, if we denote the specific 
5=1 case of the network as (it has no noniterated edges, and thus the 
maximum APL), then one can obtain F t from F^ by adding noniterated edges 
in a certain way. The less the parameter 5, the more the added noniterated 
edges. In the case of 5 = 0, the network is the densest one with the minimum 
APL. Therefore, the APL is a decreasing function of 5. As 5 drops from 1 to 
0, dt increases from ^-2* at 5 = 1 to |r2* at 5 = 0. From above arguments, 
we can conclude that for the full range of 5 between and 1, the APL dt has 
an exponential growth with network size, which indicates that the considered 
networks F t are not small worlds. 

2.5 Fractal dimension 

To determine the fractal dimension, we follow the mathematical framework 
introduced in Ref . [16] . We are concerned about three quantities, network size 
Nt, network diameter D t , and degree k u {t) of a given node u. By construction, 
we can easily see that in the infinite t limit, these quantities grow obeying the 
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following relation: N t ~ 4JV t _i, D t = 2 A-i, k u (t) ~ 2k u (t - 1). Thus, for 
large networks, Nt, k u (t) and D t increase by a factor of fx = 4, = 2 and 
fu = 2, respectively. 

From above obtained microscopic parameters demonstrating the mechanism 
for network growth, we can derive the scaling exponents: the fractal dimension 
d,B — fajP; = 2 and the degree exponent of boxes dk = yf-^ = 1. According to 
the scaling relation of fractal scale-free networks, the exponent of the degree 
distribution satisfies 7 = 1 + ^ = 3, giving the same 7 as that obtained in 
the direct calculation of the degree distribution. 

The fractality behavior can be also easy to see by tiling the system using a 
renormalization procedure as follows, we can 'zoom out' (i.e. renormalize) the 
network by replacing each connected cluster of bonds on the right of Fig. [1] 
by a single 'super'-link, in a way that reverses the process of network growth, 
see Fig. [H This has the effect of rescaling diameter, it reduces the diameter 
by a factor of 2 by carrying a cluster of bonds with diameter 2 into a single 
'super'-link with unit length. At the same time, the number nodes of in the 
rescaled network decreases by a factor = 4. Thus, the fractal dimension is 

d B = U2 = 2 - 



3 SIR model on the networks 

As demonstrated in the preceding section, the networks exhibit simultaneously 
many interesting properties: scale-free phenomenon, tunable clustering, "large- 
world" behavior, and fractality. To the best of our knowledge, these features 
have not been reported in previous network models. Hence, it is of scientific 
interest and great significance to investigate dynamic processes taking place 
upon the model to inspect the effect of network topologies on the dynamics. 
Next we will focus on studying the SIR model of epidemics, which is one of 
the most important disease models and has attracted much attention within 
physical society, see [01 and references therein. 

The standard SIR model [5|7] assumes that each individual can be in any 
of three exclusive states: susceptible (S), infected (I), or removed (R). The 
dynamics of disease transmission on a network can be described as follows. 
Each node of the network stands for an individual and each link represents 
the connection along which the individuals interact and the epidemic can be 
transmitted. At each time step, an infected node transmits infection to each of 
its neighbors independently with probability A; at the same time the infected 
node itself becomes removed with a constant probability 1, whereupon it can 
not catch the infection again, and thus will not infect its neighbors any longer. 
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The above-described SIR model for disease dynamics is equivalent to a bond 
percolation process with bond occupation probability equal to the infection 
rate A . The connected clusters of nodes in the bond percolation process 

then correspond to the groups of individuals that would be infected by a dis- 
ease outbreak starting with any individual within the cluster it belongs to. It 
was shown that in scale-free networks such as Barabasi-Albert (BA) network, 
for any vanishingly small A, there exists an extensive spanning cluster or "gi- 
ant component" , implying that scale-free networks demonstrate an absence of 
epidemic thresholds. However, we will present that for our networks, epidemic 
thresholds are present, the values of which depend on the parameter q. 

According to the mapping between the SIR dynamics and the bond percola- 
tion problem, the epidemic threshold corresponds precisely to the percolation 
threshold. In our case, the recursive construction of the networks make it pos- 
sible to study the percolation problem by using the real-space renormaliza- 
tion group technique [3"Tll3B~jl 39j to obtain an exact solution for the percolation 
(epidemic) thresholds. Here we are interested in only the percolation phase 
transition point, excluding the size of giant cluster. Let us describe the pro- 
cedure in application to the considered network. Supposing that the network 
growth stops at a time step t —>■ oo, then we spoil the network in the following 
way: for an arbitrary present link in the undamaged network, it is retained in 
the damaged network with probability A. Then we invert the transformation 
shown in Fig. 1 and define n = t — r for this inverted transformation, which is 
in fact a decimation procedure [39j. Further, we introduce the probability A n 
that if two nodes are connected in the undamaged network at r = t — n, then 
at the nth step of the decimation for the damaged network, there exists a path 
between these vertices. Here, Ao = A. According to the network construction 
and the analysis in [39] , it is easy to obtain the following recursive relation for 




Equation (jf)]) has five roots (i.e., fixed points), among which two are invalid: 
one is greater than 1, the other is less than 0. The other three fixed points are 
as follows: two stable fixed points at A = and A = 1, and an unstable fixed 
point at A c that is the percolation threshold. We omit the expression of A c as 
a function of q, because it is very lengthy. We show the dependence of A c on 
q in Fig. HJ which indicates that the threshold A c increases almost linearly as 
q increases. When q grows from to 1, A c increases from 0.5 to v ^~ 1 « 0.618. 

Therefore, there exists a percolation threshold A c such that for A > A c a giant 




A n : 




+ (l-g) A* + 5A*(1 - A„) 



+ 8A3(1-A n ) 2 + 2A^(1-A n ) 3 • 



(6) 
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Fig. 4. The dependence relation of percolation threshold A c on the parameter q. 

component appears, for A < A c there are only small clusters. This means that 
for the SIR model the epidemic prevalence undergoes a phase transition at a 
nonzero epidemic threshold A c . If the infection rate A > A c , the disease spreads 
and becomes persistent in time; otherwise, the infection dies out gradually. The 
existence of epidemic thresholds in our networks is in sharp contrast with the 
null threshold found in a wide range of stochastic scale-free networks of the 
Barabasi- Albert (BA) type [9lllUllllj . 

Why are the present scale-free networks not prone to disease propagation as 
previously studied uncorrelated BA type scale-free networks? We argue that 
the presence of finite epidemic thresholds in our networks lies with their two- 
dimensional fractal structure with diameter increasing as a square power of 
network order, a property analogous to that of two-dimensional regular lat- 
tice [10] . The 'large- world' feature stops the diffusion of diseases, and makes 
the behaviors of disease spreading in our networks similar to those of reg- 
ular lattices. Thus, the fractal topology provides protection against disease 
spreading. 



4 Conclusions 

In the present work, we have introduced a new model for fractal networks and 
provided a detailed analysis of the structural properties, which are related 
to the model parameter q. The model exhibits a rich topological behavior. 
The degree distribution is power-law with the degree exponent asymptotically 
approaching 3 for large network order. The clustering coefficient is changeable, 
which can be systematically tunable in a large range by altering q. Particularly, 
the networks are topologically fractal with a fractal dimension of 2 for all 
q. Along with the fractality, the networks display 'large-world' phenomenon, 
their average path length increases approximatively as a square power of the 
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number of nodes. 



We have also investigated the effects of the particular topological character- 
istics on the SIR model for disease spreading dynamics. Strikingly, We found 
that epidemic thresholds are recovered for all networks regardless the value of 
q, which is in contrast to the conventional wisdom that being prone to disease 
propagation is an intrinsic nature of scale-free networks. We concluded that 
the dominant factor suppressing epidemic spreading is the fractal structure ac- 
companied by a 'large- world' behavior. The peculiar structural properties and 
epidemic dynamics make our networks unique within the category of scale- 
free networks. Our study is helpful for designing real-life networked systems 
robust to epidemic outbreaks, and for better understanding of the effluences 
of structure on the propagation dynamics. 
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